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Abstract. Transport properties provide important information about the mobil- 
ity, elastic and inelastic of scattering of excitations in solids. Heat transport is well 
understood for phonons and electrons, but little is known about heat transport by 
magnetic excitations. Very recently, large and unusual magnetic heat conductiv- 

C/2 , ities were discovered in low-dimensional quantum magnets. This article summa- 

rizes experimental results for the magnetic thermal conductivity Kmag of several 

2 ' compounds which are good representations of different low-dimensional quantum 

spin models, i.e. arrangements of S=l/2 spins in the form of two-dimensional (2D) 

, square lattices and one-dimensional (ID) structures such as chains and two-leg 

^ ' ladders. Remarkable properties of Kmag have been discovered: It often dwarfs the 

usual phonon thermal conductivity and allows the identification and analysis of 
different scattering mechanisms of the relevant magnetic excitations. 



1 Introduction 
> 

[~^ , Heat transport by magnetic excitations was originally predicted in 1936 [Tj. However, it took 

lO ' almost 30 years until the first convincing experimental evidence for magnetic heat transport by 

classical spin waves was found in ferrimagnetic yttrium-iron-garnet (YIG) [2|3|4|5| . In princi- 
^^ . pie, the analysis of this magnon heat conductivity should yield valuable information about the 

f^ ' excitation and scattering of niagnons (e.g. off defects, phonons, and electrons) as is the case 

00 , for the well- understood phononic and electronic thermal conduction [B]. However, most of the 

^D ■ early experiments on YIG and following experiments on other materials [71819] were restricted 

to magnetically ordered phases at very low temperature (T < 10 K). The first signature of mag- 
netic heat transport at higher temperatures (T > 50 K) was observed for the one-dimensional 
quantum antiferromagnet KCuFs |10| . However, only the recent theoretical prediction of dissi- 
pationless heat conduction in one-dimensional antiferromagnetic Heisenberg chains |ll|12j and 
the discovery of huge magnetic contributions in the quantum spin ladder material Sri4Cu2404i 
[13|14|15] triggered intense research on the heat transport of low-dimensional quantum spin sys- 
tems [llp :2 13 14 15 16 17 18 19 20 21 22 23 24 2 5I26|27|28|29| 30|31 32 33 34 35 36 37 38 39 40] 
[41I42I43I44.45.46.47.48.49.50.51.52.53.54 55 . Over the course of this work, more and more 
cases for low dimensional magnetic heat conduction were observed in various materials. Today, 
the clearest experimental examples of low dimensional magnetic heat conduction are found in 
copper oxide (cuprate) systems. The overview of the experimental research status on low dimen- 
sional magnetic heat conduction provided in this article therefore focuses on these compounds. 
Particular examples from the plethora of possible spin structures in cuprate systems are 
spin arrangements in the geometrical form of chains, so-called two-leg ladders, and square 
lattices with a strong antiferromagnetic Heisenberg exchange (J « 1500 — 2000 K) between 
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Fig. 1. Illustration of low-dimensional spin structures: (a) a spin chain, (b) a two-leg spin ladder, and 
(c) a two-dimensional square lattice. Arrows represent localized S = 1/2 spin and shaded bars symbolize 
strong antiferromagnetic exchange between them, (d) Schematic illustration of the underlying chemical 
building block giving rise to the spins and their interaction. Only the relevant Cu 3d^2_y2 and O 2p^ 
orbitals are indicated. Arrows represent the spins of the electrons involved. 



nearest neighbor spins. Sketches of such spin arrangements are shown in Fig.[T^-c. These low- 
dimensional spin structures usually arise from similar low dimensional Cu-0 structures, in 
which the antiferromagnetic exchange originates from straight Cu-O-Cu bonds (bonding angle: 
180°) as depicted in Fig[T]i. All these systems are based on Cu^^-ions and therefore represent 
S = 1/2 systems with a strong quantum nature. 

Good examples for materials containing S = 1/2 Heisenberg chains as depicted in Fig. [T^ 
are given by the compounds CaCu203, SrCu02 and Sr2Cu03, where straight Cu-O-Cu bonds 
and hence a strong antiferromagnetic exchange only exist along one particular crystallographic 
direction; the magnetic exchange perpendicular to this direction is much weaker |56t57j . In the 
(Sr, Ca, La)i4Cu2404i family of compounds, parallel pairs of such chains are coupled to each 
other via bridging 0-ions, producing in straight Cu-O-Cu bonds perpendicular to the chain 
direction. The resulting magnetic interchain coupling perpendicular to the chain direction J± 
is of a similar magnitude to the intrachain coupling, i.e. J± ^ J [58j . This spin ensemble is a 
so-called two-leg spin ladder, the ladder legs being formed by the two chains and the ladder 
rungs arising from the Cu-O-Cu bonds which connect the chains (cf. Fig. [T]]). Following this 
concept, ladder structures with more legs can in principle be created by coupling more chains to 
the structure; eventually this would lead to a two dimensional Heisenberg antiferromagnet on a 
square lattice (2D-IIAF) in the infinite limit. A good realization of a 2D-HAF with S — 1/2 is 
given by La2Cu04 and other related antiferromagnetic parent compounds of high-temperature 
superconductors. 

The corresponding low-dimensional quantum spin models are characterised by very peculiar 
ground state properties and elementary excitations, which vary strongly from system to system. 
Spin-spin correlations of homogeneous spin chains, for example, are quasi-long range in the 
ground state and decay algebraically with distance between the spins [59]. The elementary 
excitations, so-called spinous, are gapless and carry a spin 5* = 1/2 [60]. In contrast to this, 
the ground state of a two-leg ladder is a spin liquid, i.e. the spin-spin correlations are short 
range and decay exponentially as a function of distance. The elementary excitations are 5* = 1 
particles (usually called magnons or triplons) and possess a spin gap A {A/ks ~ 400 K in 
the case of the systems discussed here) [S^ • Finally, the ground state of the 2D-IIAF is a Nccl 
state with strongly reduced sublattice magnetization, which only exists at temperature T = 
[61j . In this case the elementary excitations are well described using a spin wave framework. 
However, it should be noted that alternative descriptions have been discussed |62|63|64j . 

In the case of hole doping, all these model systems yield interesting and exotic properties. A 
Luttinger liquid forms in hole-doped spin chains, i.e. electronic excitations decay into collective 
excitations of holes (holons) and spins (spinous) . This phenomenon is usually called spin-charge 
separation. Radically different properties have been predicted for two-leg spin ladders: super- 
conductivity competing with a charge ordered ground state is expected in this case I65I66I . 
Finally, hole doping has great importance in the case of the 2D-HAF as is evident from the 
high temperature superconductivity which is observed in such systems. Note that controlled 
hole doping of S* = 1/2 chains and the observation of spin-charge separation signatures has not 
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Fig. 2. Anisotropic thermal conductivity of various low-dimensional spin materials as a function of 
temperature: (a) the spin chain compound CaCu203, (b) the two-leg spin ladder material Sri4Cu2404i, 
and (c) the 2D-HAF as realized in La2Cu04. Filled and open symbols represent Ky and k±_ of the 
materials. The solid line in (a) represents a linear fit to the data in the range T > 100 K. The axis 
intercept of its extrapolation towards T = is an approximation of Kph in the fit range. Solid lines in 
(b) and (c) represent estimations for the phonon background. From Refs. [23115117] . 



yet been achieved experimentally, whereas charge ordering and superconductivity are prominent 
experimental features of hole-doped spin ladder and 2D-HAF materials. 

Concerning heat transport, little is known for all these model systems. Often the attention 
in theoretical works is focussed on the possibility of ballistic magnetic heat transport in ID- 
systems: in integrable models like the XXZ Heisenberg spin chain the Hamiltonian and the 
thermal current operator commute, i.e. once a thermal current is established in such a system it 
will never decay [11| . In other words, the thermal resistance vanishes and the magnetic thermal 
conductivity Kmag diverges. While such surprising properties are well established for integrable 
spin models [11], ballistic heat transport in non- integrable quasi ID-systems (e.g. two- leg spin 
ladders) is currently a subject of intense discussion ^38 46 39 55]. However, in real materials 
scattering processes involving defects and other quasiparticles such as phonons and charge 
carriers must play an important role and render Kmag finite in all cases [43] . The analysis of 
Kmag should hence provide further insight into the nature of these scattering processes and the 
dissipation of magnetic heat currents. 



2 Experimental signatures of magnetic heat conduction 



Fig. [2] shows experimental results [14115117123] for the thermal conductivity k of CaCu203, 
Sri4Cu2404i and La2Cu04 which are good experimental representations of 5' = 1/2 isotropic 
antiferromagnetic Heisenberg spin chains, two-leg spin ladders and the 2D-HAF, respectively. 
The experimental thermal transport properties of all these electrically insulating^ are materials 
are intriguing: when k is measured perpendicular to the low dimensional structure (k_l), i.e. 
a direction along which the magnetic coupling is negligible, the T-dependence of ordinary 
phonon thermal conductivity ^ Kph is found for all three materials: kj_ exhibits a peak at 
low temperature T w 20 K, which is followed by a continuous decrease as T rises further. 
Note the exception in CaCu203, where one component of k,±_ increases monotonically with 
rising T. Here, a strong suppression of Kph due to disorder and possible contributions from 
optical phonons could give rise to the observed T-dependence [23] . The situation is completely 
different for k parallel to the low-dimensional structures, i.e. along the directions with large J 
(kii). Again, a phononic low-T peak is observed. However, nn evolves very differently at higher 
T. At T > 75 K, K|| strongly increases upon heating and exhibits a peak for Sri4Cu2404i and 
La2Cu04 at 140 K and 310 K, respectively, while for CaCu203 the increase continues up to 



^ Despite the spin ladder compound Sri4_xCaxCu2404i being intrinsically hole doped, electronic 
contributions to k, are negligible. 
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Fig. 3. Estimated magnetic thermal conductivity as a function of temperature of (a) the spin chain 
compound CaCu203, (b) the two-leg spin ladder compounds Sri4Cu2404i (•) and Ca9La5Cu2404i 
(o), and (c) the 2D-HAF La2Cu04. The solid lines represent fits according to the respective expres- 
sions for Kmag in selected T-ranges (see text) with the approximate T-dependence indicated. From 
Refs. |23I15I17| . 

the highest temperature measured. In these three cases, the remarkable anisotropy of k is the 
quahtative evidence for large magnetic contributions to Ky, i.e. for magnetic heat conductivity 
Kmag- A similar anisotropy is also present for ku of the spin chain materials (Sr, Ca)Cu02 
and Sr2Cu03 which led to the conclusion that spinon heat transport is also present in these 
materials i25|26 21J . However, as will be discussed in SectionlU a high-T peak is absent in those 
cases which renders a quantitative analysis of phononic and magnetic contributions to kii more 
difficult. However, a clear separation of kii into magnetic and phononic parts (Kmag and Kph) is 
possible for the three cases shown in Fig. [2] where the strong features of Kmag appear at a much 
higher T-scale than the low-T phonon peak. In the following sections we will therefore focus on 
these cases and examine what can be learned. 



2.1 Extraction of magnetic contributions 

In order to obtain the magnitude and T-dependence of Kmag, it is essential to accurately estimate 
Kph and subtract it from the total K||. In the case of La2Cu04 and Sri4Cu2404i this can be 
performed in a convenient way since the magnetic contributions are expected to be negligibly 
small in comparison to Kph for T < 40 K, i.e. in the vicinity of the phononic peak. This is 
due to the expected T-dependence of Kmag being approximately ex T^ and oc exp —A/T for 
the 2D-HAF and the spin ladders, respectively (cf. also the analysis further below). We may 
hence estimate K|| ~ Kph at T < 40 K and to extrapolate the low temperature Kph to higher 
T. The thus estimated Kph is indicated in Fig. ^Bp and [2J: as solid lines. For details of the 
procedure the reader is referred to Refs. [14I15I17J . In the case of CaCu203 this procedure is 
not applicable since Kmag oc T is expected for a 5 = 1/2 Heisenberg chain hence providing 
significant contributions to kii even at low T. However, in the present case of CaCu203 the 
situation is quite fortunate, since K|| 3> k_l ~ /^ph ~ const at T > 100 K which allows the 
extraction of Kmag from the total K|| at temperatures higher than 100 K simply by subtracting 
a constant value j23]. For all three cases the conjectured Kmag are shown in Fig. [S] The figure 



also shows 



•rnag 



of La5Ca9Cu2404i, which belongs to the same family of two-leg spin ladder 



compounds as Sri4Cu2404i but differs by a lower content of charge carriers in the ladders (cf. 
discussion in Section [3731) . 



3 Analysis of magnetic heat conductivity 

We start our analysis by considering the qualitative T-dependence of Kmag, which comprises a 
simple peak structure (La2Cu04 and Sri4Cu2404i) and a monotonic increase (CaCu203) in the 
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studied range T — 100-350 K, where the latter may be regarded as the low temperature edge 
of a peak. A peak structure is very common for the thermal conductivity k of any kind of heat 
carrying particle, such as phonons or electrons 6 . In the following we will investigate whether 
the underlying physics can be applied to magnetic excitations as well. The basic physics which 
determines the T-dependence of k can be inferred from the kinetic estimate 

' CkWk^krfk, (1) 



d (27r)" 

with d the dimensionality of the considered system, Ck — ^eic'^^k the specific heat (ek and rik 
are the energy and the statistical occupation function of the mode k), vit_ the velocity and ?k the 
mean free path of a particle with wave vector k. At low T only a few particles are excited and 
contribute to the heat transport. Often scattering processes are rare and l^ is a slowly varying 
function of momentum in the relevant energy range. In this situation, the low-T increase of k is 
(a) characteristic of the excitation of the heat carrying particle (reflecting the T-dependence of 
the specific heat if Uk is momentum independent) and (b) proportional to the mean free path 
I ~ Ik. At higher T the momentum-dependent scattering becomes more important and leads 
to a decrease of the mean free path and hence to a decrease of k. Normally this decrease is 
characteristic of the relevant scattering mechanisms and allows an advanced analysis. 

The application of Eq. [1] for the case of ID and 2D magnetic systems leads to the general 
result 

Kmag(T) OC l„,agf{T) , (2) 

where lma.g{T) is a general magnetic mean free path based on the approximation Zmag = ^k- For 
fc^T ^ J this assumption is justified because the heat carrying excitations exist in significant 
numbers only in the vicinity of the band minima, i.e. a very small fraction of the Brillouin zone. 
This requirement is always fulfilled for all systems discussed in this article since J « 1500- 
2000 K and the experimental data only extend over temperatures T < 350 K. The function 
/(T) depends on temperature in a manner which is characteristic of the considered spin systemo 
In particular, for a gapless S = 1/2 Heisenberg chain Eq.[T] yields [55] 



where Hs is a geometrical factor that counts the number of chains per unit area. At low tempera- 
tures ksT <^ J the integral is only weakly temperature dependent and approaches the constant 
value 7r^/6 for T — > 0. Note that the condition kgT ^ J holds even at room temperature (up 
to T < 0.15J//cB « 300 K), i.e. the experimental data always represent the low temperature 
behaviour of Kmag [26 23] and the upper boundary of the integral may be set to infinity. The 
corresponding expression for a gapped two-leg ladder is [15j 

_ 3n,kB^ f°° 2 exp(x) 

«mag - ^^ imagJ J ^_ X ^^^^^^^ ^ ^^^ax , [i) 

kgT 

where Us is now the number of ladders per unit area. At low temperatures ksT <^ A < J one 
might approximate /(T) ex exp{—A/{kBT)). Finally, for a 2D-IIAF one finds, accounting for 
both magnon branches [17j . 

_ Y^ TipkB^ , „2 /"°° 2 / 2 2 exp(x) ^ 

Km.s - 2^ 4^;^2„p '"-S 7^ ^ V ^ ""O-' (exp(x) - 1)^ """^ ' ^^' 

i — 1,2 k^ T 

with the spin wave velocity vq « 1.287 • 10^ m/s ^7\ and Up = 2/c the number of planes 
per unit length along the c-axis, which is the direction perpendicular to the planes (c is the 



^ For all three different types of systems the reader is referred to the original literature [15117126123) 
for the derivation of the respective f{T). Within this article, we specify only the final results for Kmag- 
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corresponding lattice constant). The integral is temperature dependent via its lower boundary 
xo.i — Ai/{kBT), where Ai/ks « 26 K and A2/kB ~ 58 K account for the anisotropy gaps 
which arise in La2Cu04 [68| . However, the temperature dependence is weak in the T-range 
where the experimental data are discussed and one can approximate K„iag oc T^ . 



3.1 Low-temperature characteristics - thermal occupation 

Interestingly, the experimental data shown in Fig. [3] exhibit extended regions at low T where 
a reasonable description of the data with Eq. [2] using a temperature independent mean free 
path Zmag is possible. The solid lines in Fig. [3] represent fits where /mag is a free, temperature 
independent parameter. A remarkably good description is found for Kmag of the spin chain 
CaCu2 03, as depicted in Fig. [3^. Here, n-axag is excellently described by a simple linear increase 
over the large temperature range T k, 100-300 K and the fit yields Zmag ~ 22±5 A corresponding 
to about 5-6 lattice spacings. 

For the spin ladder compounds (Sr, Ca, La)i4Cu2404i the situation is somewhat different. 
As can be seen in Fig. [3}d, the temperature range over which a good fit with Eq. |4] can be 
achieved is strongly reduced with respect to the previous case and apparently also depends on 
the composition of the material. In particular, the temperature interval where Eq. 3] describes 
K,„ag with a T-independent mean free path is 54-102 K for Ca9La5Cu2404i but only 61-91 K 
for Sri4Cu2404i. Nevertheless, restricted to these T-ranges the fit (with the spin gap A and 
the mean free path Z^ag as free fit parameters) yields similar results for both compounds. 
In particular, the value for the spin gap is found to be somewhat larger than but still in 
reasonable agreement with spin gap results from neutron scattering {A/kB = 418 ± 15 K 
and A/kB = 396 ± 10 K for Ca9La5Cu2404i and Sri4Cu2404i, respectively) J69I70I71| . The 
magnetic mean free path is of similar magnitude in both materials (/,nag = 2980 ± 110 A and 
^mag = 2890 ± 230 A respectively) O It is important to note that for temperatures higher than 
the mentioned ranges, Eq. |4] completely fails to properly describe the experimental data with a 
T-independent ^mag, as is also evident from Fig.[3)D. However, as will be discussed in Section [373l 
a consistent picture arises if ^mag is allowed to become T-dependent at higher T. 

A similar observation is also made in the case of the 2D-HAF La2Cu04 (cf. Fig. [3t). Eq. [5] 
with Zniag as a free parameter yields a reasonable fit in the range 70-158 K, where /mag ~ 
558 ± 140 A [ITj. Again, the theoretical model fails to account for the high temperature regime 
as long as /mag remains T-independent. 



3.2 Low temperature characteristics - scattering off defects 

The actual meaning of the magnetic mean free path /mag as a material parameter is a priori not 
clear. On one hand it is known from the analogous case of phonon heat transport, that at low 
temperature the phonon mean free path can become as large as the crystal dimensions, i.e. of the 
order of millimeters [72]. However, the extracted values for /mag are several orders of magnitude 
smaller than the dimensions of the crystals which have been studied in the experiments. It 
therefore appears natural to conclude that /mag should reflect the density of static defects in 
the material. On the other hand, a magnetic mean free path of the order of up to ^ 1000 lattice 
spacings implies almost perfect crystallinity of the underlying material, which is astonishing 
in view of the large and complicated unit cell of, for example, Ca9La5Cu2404i. Note that a 
much smaller general magnitude of /mag has been suggested for this compound based on Exact 
Diagonalization calculations of the thermal Drude weight which implied ballistic heat transport 
in spin ladder systems |38| . Such reduced values for /mag could not be confirmed by more recent 
calculations which suggest a vanishing Drude weight [39146) . 



^ The slightly smaller values of /mag and A for Ca9La5Cu24 04i as compared to previous results [15] 
are a consequence of the usage of more accurate lattice parameters and an optimized fit-interval. It is 
stressed that these small corrections have no further consequences on the conclusions drawn in Ref. [15] . 
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Fig. 4. a) Magnetization of CaCu203 as a function of temperature with a magnetic field [iqH — 1 Tesla 
parallel to the three crystallographic axes. The solid lines represent Curie- Weiss-type fits to the data 
in the range 100-360 K which yield ~ 3% free spins with respect to Cu. Inset: Inverse magnetization 
after subtracting a constant A/q which roughly accounts for van Vleck and chain magnetism. From 



Ref. |23) . b) Main panel: Magnon thermal conductivity 



of La2Cui_^Zn^04 {z = 0, 0.005, 0.008, 



0.01, 0.02, 0.05) as a function of T (open circles). Solid lines: fits according to Eq. [5l Inset: /mag as a 
function of 1/z in units of lattice constants a. Solid line: fit through origin. From Ref. 17 . c) /mag of 
Sri4Cu24-zZnz04i (z = 0.125, 0.25, 0.5, 0.75) as a function of the mean distance between Zn impurities 
dzn-Zn- Solid line: linear fit through the origin. From Refs. |23I22I17] . 



A straightforward experimental method to elucidate the connection between ^^ag and the 
density of magnetic defects in the material is to measure the latter independently from heat 
transport using a different experimental technique. One possibility is to study the magnetic 
susceptibility x(^) where one has to assume that paramagnetic moments, the concentration of 
which is deducible from x, are connected with defects within the magnetic structure and hence 
may have an effect on the magnetic heat transport. For CaCu203 the situation appears to be 
quite fortunate, since in a recent study Goiran et al. suggested a direct link between paramag- 
netic moments located off the ID magnetic chain structures (detectable by x-measurements) 
and (possibly non-magnetic) defects within the chains [73]. From susceptibility measurements 
on CaCu2 03 (cf. Fig. [4]) it was then possible to deduce an upper limit for the mean distance 
between defects within a chain which turned out to be a factor of 2-3 larger than the extracted 
^mag [23j . The same order of magnitude of both quantities indicates that these defects are the 
main scatterers for magnetic excitations within a chain. In the absence of a similar method 
for La2Cu04 and (Sr, Ca, La)i4Cu2404i, an alternative approach was selected to compare /mag 
with known distances between intentionally doped defects in the material. Such defects can be 
induced by substituting a small amount of non-magnetic Zn^+ for the magnetic Cu^^-ions. This 
has been performed for both La2Cu04 and Sri4Cu2404i for a series of different doping levels. 
As can be inferred from FiglDa andHJ:, once again /mag has the same order of magnitude as the 
mean distance between the defects, which is in these cases well defined by the mean distance 
between the Zn dopants [17I22J . This quantitatively confirms that a good analysis of Kmag in 
these low dimensional spin systems can be performed using a rather simple kinetic model. 



3.3 Scattering processes at higher temperatures 



We now turn briefly to the behavior of Kmag at high T in order to elucidate the impact of 
temperature dependent scattering processes on the magnetic heat transport, i.e. scattering of 
magnetic excitations off other quasiparticles such as phonons, charge carriers or other magnetic 
excitations. Here we focus on heat transport in the spin ladder compounds, since these are 
susceptible to a variety of different types of doping. Most remarkable (apart from doping with 
non-magnetic impurities which has already been discussed) is certainly the possibility of hole- 
doping the ladders. Interestingly, the stoichiometric parent compound Sri4Cu2404i is already 
inherently doped with holes. These holes reside only partially in the two- leg ladder structures 
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- the largest portion is located in chain substructures which are also present in the material 
beside the two-leg ladders. The holes are redistributed between the chain and ladder structures 
upon the isovalent substitution of Ca for Sr: with increasing Ca-content a significant increase 
of the hole concentration in the ladders is observed |74|75|76] . It is also possible to reduce 
the overall hole content in the material and thereby render the ladders virtually hole-free by 
replacing the divalent Sr or Ca by trivalent ions such as La [74I75J . 

A good example for the latter case of doping is the compound CagLa5Cu2404i whose mag- 
netic thermal conductivity is shown in Fig. [8)3 in comparison with that of Sri4Cu2404i |18j . 
The effect of hole-doping on Kmag can immediately be observed in this figure. For T < 100 K 
the increase in K^ag with T is almost identical for both compounds. Pronounced differences 
only occur at higher T: Kmag of La5Ca9Cu2404i exhibits a large peak (~ 140 Wm^^K^^ at 
~ 180 K) and stays very large even at room temperature (~ 100 Wm^^K^^). In contrast, the 
peak is much smaller in the case of Sri4Cu2404i (~ 75 Wm~^K~^ at ~ 150 K). Here K^ag 
decreases much more strongly at high T and saturates at Kmag ~ 10 Wm^^K^^ for T > 240 K. 

It is straightforward to attribute the strong high-T suppression of Kmag in Sri4Cu2404i 
compared with La5Ca9Cu2404i to scattering of the magnons off holes, since the hole dop- 
ing in this compound is the most relevant difference with respect to the undoped ladders 
of La5Ca9Cu2404i. Both K^ag curves are almost identical below a characteristic temperature 
To « 100 K, confirming that this scattering mechanism becomes completely unimportant below 
To and only reveals its full strength above a characteristic temperature T* « 240 K. 

The surprising temperature dependence of the strength of the magnon-hole scattering has 
been checked for robustness against changes of the hole content in the ladders, where Kmag of a 
series of Sri4_xCaxCu2404i single crystals has been investigated in detail [18]. The comparison 
with Kniag of the Ca-doped samples (see Fig. [5^) reveals that Tq and T* are gradually shifted 
towards lower T; i.e., the temperature region where K^ag is suppressed extends and magnon- 
hole-scattering also becomes important at low T. At a; = 4, 5 this region appears to becomes so 
wide that even the peak at low T is suppressed. It was shown that the temperature dependence 
of Kniag is unambiguously correlated with a charge ordered state in the compound, where charge 
ordering sets in at T < T*. More precisely, the charge ordering in the ladders is accompanied 
by a drastic enhancement of the magnon mean free path Imag- the probability for magnon-hole 
scattering, which is close to unity for mobile holes, vanishes in the charge ordered state [18] . 

Turning to Kmag of La5Ca9Cu2404i, i.e. that of undoped ladders, it is a priori clear that the 
magnon hole scattering, which is dominant in Sri4_xCaxCu2404i, cannot play any significant 
role. Nevertheless, the necessity to allow a formal T-dependence of /,nag at T > 100 K indicates 
that further scattering processes are also relevant for this case. The only remaining possible 
scattering mechanism for magnons in this material are magnon-magnon or magnon-phonon 
scattering. In a more involved analysis the T-dependence of /mag in La5Ca9Cu2404i has been 
calculated from the K^ag data using Eq. [3]by employing the previously extracted A = 418 K. 
The resulting /mag(T) as shown in the main panel of Fig. [Sb reflects the different T-regimes 
which govern Kmag- For T < 110 K, Imug is T-independent with a mean value Iq — 2980 A which 
reflects the scattering of magnons off static defects. In order to describe the T-dependent Zmag 
at higher T it was assumed that all scattering mechanisms were independent of each other and 
Matthiesen's rule applied: Z~ag = ^o" + 7phrfph + 7magrfmig- Here dph and dmag are the mean 
" distances" of phonons and magnons respectively, as calculated from the particle densities with 
7ph and 7mag the corresponding scattering probabilities. Since it is unclear as to what extent 
the separate scattering mechanisms contribute to Imag, its behavior was analyzed based on the 
assumption that only one mechanism is active in addition to magnon-defect scattering. 

The case of dominant magnon-phonon scattering was modelled by three energy-degenerate 
non-dispersive optical branches along the ladder direction, yielding 1/dph = 7.6 ,^ ^ /t)-i ' 

10^m~^ with Z\opt the optical gap (cf. Ref.[20] for details). The experimental Imag, was then 
fitted with /~ag — ^a^ + Tph^^ph using 7ph and Z\opt as free parameters. The fit (solid line in 
Fig. [Sb) describes the data fairly well. Remarkably, the value found for Z\opt = 795 K is of 
the same order of magnitude as the energy of the longitudinal Cu-0 stretching mode which is 
involved in the two-magnon-plus-phonon absorption observed in optical spectroscopy [771 78j . 
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in Fig.tSb). cl is the lattice constant along the ladders and e^ was taken from Johnston et al. for 
the case of isotropic ladder coupling 79J, with £k=TT = A — 418 K employed. The comparison 
between both fits suggests that scattering off optical phonons is dominant in this compound. 



4 Other developments 

Complementary to Kmag of the "dirty" spin chains of CaCu2 03, the spinon heat transport of 
the very "clean" spin chain materials Sr2Cu03 and SrCu02 has been investigated by A. V. 
Sologubenko and coworkers [25|26j . Despite the data giving very clear evidence for spinon heat 
conductivity in this compounds, the precise extraction of Kmag from the experimental data is 
quite involved, since the signature of Kmag appears as a shoulder in the high-T edge of the 
phononic low-T peak (see Fig. [5fc for the case of SrCu02). In their analysis, Sologubenko ct al. 
suggest that Kmag exhibits a peak-like T-dependence which allows for further analysis 25 26J. 
In order to achieve a better extraction of Kmag, P- Ribeiro et al. investigated the thermal 
conductivity of Sri_xCaxCu02 with x — 0.05 with the expectation that the Ca-impurities 
selectively suppress Kph while Kmag remains unchanged [H]. However, as can be seen in Fig.jnt, 
the Ca-doping apparently leads to a complete suppression of the shoulder-like anomaly of the 
phonon peak and no signature of a peak structure remains (cf. Ref. [21 for details). Further 
studies are necessary to elucidate the origin of this intriguing observation. 



5 Conclusion 



From this summary of recent developments in the research on magnetic heat transport in quan- 
tum spin systems, it becomes evident that even though considerable progress has been made, 
this new thermal transport mechanism is still far from being thoroughly understood. We have 
seen that in some cases the magnetic heat conductivity can serve as a sensitive probe for mag- 
netic excitations, which is a promising new approach to access the scattering mechanisms and 
dissipation of these excitations. However, many issues still remain to be resolved. For example, 
little is known about how Kmag evolves when the magnetic systems become less quantum in 
nature, i.e. when S > 1/2. Initial experiments and theoretical work have already addressed 
this topic |33|47|80(54] but the number of investigated materials of this type is still small. A 
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better understanding of scattering processes (non- magnetic vs. magnetic impurities, phonons), 
frustration, and the effect of an external magnetic field will also be required. 

This overview was only made possible by the valuable contributions of P. Ribeiro, C. Baumann, H. 
ElHaes, A. Waske, U. Ammerahl, L. Colonescu, G. Krabbes, P. Rentier, A. Revcolevschi, G. Roth, C. 
Sekar, F. Heidrich-Meisner, W.Brenig, and B. Biichner. It is a pleasure to further thank A. Chernyshev, 
J. Gemmer, A. Honecker, M. Michel, A. Rosch, A.V. Sologubenko, and X. Zotos for fruitful discus- 
sions, and A. P. Petrovic for proofreading the manuscript. This work was supported by the Deutsche 
Forschungsgemeinschaft through grant HE3439/7 and the European Commission through the FET 
open-STREP NOVMAG, Project Reference 032980. 
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